Matrix Representation of Renormalization in 
Perturbative Quantum Field Theory 



KURUSCH EBRAHIMI-FARD 
Physics Institute, Bonn University, Nussallee 12, 
Bonn 53115, Germany 

e-mail: f ardSth . physik . uni-bonn . de 

and 
LI GUO 

Department of Mathematics and Computer Science, Rutgers University, 

Newark, NJ 07102, USA 

e-mail: liguo@newark.rutgers.edu 

August 20, 2005 



Abstract 

We formulate the Hopf algebraic approach of Connes and Kreimer to renormalization in 
perturbative quantum field theory using triangular matrix representation. We give a Rota- 
Baxter anti-homomorphism from general regularized functionals on the Feynman graph Hopf 
algebra to triangular matrices with entries in a Rota-Baxter algebra. For characters mapping 
to the group of unipotent triangular matrices we derive the algebraic Birkhoff decomposition 
for matrices using Spitzer's identity. This simple matrix factorization is applied to characterize 
and calculate perturbative renormalization. 
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1 Introduction 



Most, if not all, of the interesting and relevant 4-dimensional quantum field theories suffer from 
ultraviolet divergencies and need to be renormalized [f5]. The basic idea of the theory of perturba- 
tive renormalization in quantum field theory goes back to Kramer [12], and was successfully applied 
for the first time in a 1947 paper by Bethe [6], dealing with a concrete problem in perturbative 
quantum electrodynamics (QED). Five decades later. Dirk Kreimer [35] uncovered a Hopf algebra 
structure underlying the intricate combinatorial-algebraic structure of renormalization in general 
perturbative quantum field theory (QFT), hereby providing a sound and useful mathematical foun- 
dation for this most important achievement of theoretical physics. Later, Kreimer [36, 37, 38, 39] 
and collaborators [5, 7, 8, 9, 10, 11, 22], especially Connes and Kreimer [16, 17, 18, 19] further 
developed the Hopf-algebraic approach, connecting it to non-commutative geometry. 

In their approach, one particle irreducible (IPI) Feynman graphs, as the building blocks of 
perturbative QFT and renormalization, are organized into a connected, graded, commutative, non- 
cocommutative Hopf C-algebra T-Lj^. The restricted dual of this Hopf algebra of Feynman graphs, 
denoted by TV-p, contains the group G :— char{Ti.j^, C) of characters, that is, algebra homomorphisms 
from TijF to the underlying base field C. Feynman rules naturally provide a special class of such 
characters. The group G is generated by the Lie algebra g :— dchar{Hj^, C), formed by derivations, 
or so-called infinitesimal characters. We refer the reader to [26, 29, 37, 40, 49] for more details. 

Dealing with the ultraviolet divergencies demands a regularization plus a renormalization scheme. 
As a main example for the former serves dimensional regularization, where we replace the above 
base field C by Laurent series A — C[e~^,e]\. They form equipped with the ordinary multipli- 
cation a commutative Rota-Baxter algebra, and we consider Ga '■— char{Hj^, A), respectively 
qa '■— dchar{Hj^,A). The Hopf algebra of Feynman graphs and, via the famous Milnor-Moore 
Theorem [29, 43], its equivalent Lie algebra and Lie group of characters, enabled Connes-Kreimer 
to capture the process of renormalization in terms of a so-called algebraic Birkhoff decomposition of 
regularized Feynman rules characters, giving rise to a link to the Riemann-Hilbert problem [17, 18]. 
See also [20]. 

In [24, 25] the intimate link between the notion of Rota-Baxter algebras and the work of Connes- 
Kreimer was explored. It was shown that a non-commutative generalization of a classical result 
with origin in fluctuation theory of probability, known as Spitzer's identity for Rota-Baxter algebras 
[3, 47], lies at the heart of the Connes-Kreimer decomposition theorem for regularized Hopf algebra 
characters. This provides a natural way to derive Bogoliubov's recursions for the counter term 
and renormalized Feynman rules. This approach emphasizes the Lie algebra structure on Feynman 
graphs and the corresponding Lie group of characters. The former is closely related to the more 
general insertion and elimination Lie algebra of Feynman graphs, which was studied in detail first 
by Connes-Kreimer in [19], and recently in the context of matrices by Mencattini and Kreimer 
[41, 42]. 
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The Hopf algebras of Feynman graphs found in renormahzation theory fall into the class of so- 
called combinatorial Hopf algebras. In this paper, we obtain a matrix representation of the full space 
Hom(7i^^, A) of regularized functionals, thus in particular of Ga and qa, directly from the coproduct 
of the Hopf algebras Hj^. In fact, the matrix representation applies to a large class of Hopf algebras H 
and gives an injective anti-homomorphism of Rota-Baxter algebras from Hom(7i, A) to the infinite 
size upper triangular matrices Ai'^{A). Therefore, in the particular case of QFT renormahzation, 
i.e., when 7i — Tijr, the whole process of perturbative renormahzation of regularized Feynman 
characters, presented in the context of Rota-Baxter algebras [24, 25], is translated to a parallel, but 
more transparent process involving matrices. The upper (or lower) triangular matrices with entries 
from the target space of such Hopf algebra functionals form a complete filtered non-commutative 
associative Rota-Baxter algebra, giving rise to Spitzer's identity and a matrix factorization, from 
which wc easily read off the Birkhoff decomposition of Connes-Krcimcr. 

Our work was partly motivated by Berg and Cartier's [4], which used the pre-Lie insertion 
product on Feynman graphs to investigate the renormahzation Hopf algebra of Connes and Kreimer 
in terms of a lower triangular matrix representation of the Lie group Ga and its generating Lie 
algebra qa- We will compare with their work in subsection 3.6. 

The following is a summary of the paper. In the next section we review the notion of complete 
filtered Rota-Baxter algebras and Spitzer's identity. Theorem 5 is a generalization of Spitzer's 
classical result to associative non-commutative Rota-Baxter algebras, giving rise to a factorization 
theorem with recursively defined solutions. As an interesting and useful application of this theorem 
we formulate explicitly the factorization of unipotent triangular matrices with entries in a com- 
mutative Rota-Baxter algebra, and provide detailed calculations of small size matrices. Section 3 
deals with the application of the matrix Rota-Baxter algebras to the Hopf-algebraic approach to 
the process of renormahzation in perturbative QFT. Subsection 3.2 is the heart to this paper. We 
define a representation of the dual space of the renormahzation Hopf algebra of Connes-Kreimer in 
terms of upper triangular matrices, using the coproduct structure map. This gives rise to a (Lie) 
group of upper triangular matrices with unit diagonal, representing the group of regularized Hopf 
algebra characters. The associated generating Lie algebra of infinitesimal characters, or derivations, 
is represented by nilpotent upper triangular matrices. We include several explicit examples, calcu- 
lating the renormahzation of amplitudes of Feynman diagrams up to three loops in dimensionally 
regularized four dimensional (/7^-theory as a quantum field theory toy model. Further applications 
and calculations of our results to renormahzation in perturbative QFT are detailed in the compan- 
ion article with J. M. Gracia-Bondia and J. C. Varilly [23]. We finish with a comment on the work 
of Berg and Cartier giving matrix Lie algebra representation using the pre-Lie algebra structure on 
Feynman graphs. 
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2 Complete Rota— Baxter algebras of triangular matrices 

In this paper K denotes a field of cliaracteristic zero witli unit denoted by Ir and often by 1. All 
algebras are assumed to be unital associative K-algebras, if not stated otherwise, with unit identified 
with Ik- 

2.1 Rota— Baxter algebras 

We recall basic concepts and properties of Rota-Baxter algebras. The reader may consult the 
following literature [2, 3, 30, 31, 44, 45, 46] for more information. 

By a Rota— Baxter algebra in this paper, we always mean a Rota-Baxter K-algebra of weight 
one, that is, a K-algebra A with a Rota-BcLxter map R : A, fulfilling the relation 

R{x)R{y) + R{xy) = R{R{x)y + xR{y)) , Vx, e A. (1) 

So we also denote a Rota-Baxter algebra by a pair [A,R). For later reference we mention the 
example of the well-known dimensional regularization scheme where A = C[e~^,e]], the field of 
Laurent series. One easily shows that A is a commutative Rota-Baxter algebra, with the Rota- 
Baxter map R := R^s being the pole part projection, known under the name of minimal subtraction 
(MS) scheme 

Rms I ^ I '■— ^k^^- 

\k=-N / k=-N 

For a Rota Baxter map R on A, the map R :— id^ — i? is also a Rota-Baxter map. Further we 
have the mixed relation 

R{x)R{y)=R{xR{y))+R{R{x)y), Vx,y e A. (2) 

The images of both Rota-Baxter maps R and R form non-unital subalgebras of A. As a trivial 
observation we state the fact that every algebra A is a Rota-Baxter algebra with Rota-Baxter pair 
id^ and id^ = 0. A homomorphism / : {A, R) {A', R!) between two Rota-Baxter algebras is a 
ring homomorphism such that f o R = R' o f. 

A Rota— BcLxter ideal of a Rota-Baxter algebra {A, R) is an ideal I oi A such that R{I) C /. 
Let {A, R) be an associative Rota-Baxter algebra. The Rota-Baxter relation extends to the Lie 
algebra Ca with commutator [x, y] :— xy — yx, Vx, y e A. In other words, 

[Rix),Riy)] + R{[x,y]) = R{[R{x),y] + [x,R{y)]),y x,y e Ca: 

making {Ca, R) into a Rota— Baxter Lie algebra. 

Given a Rota-Baxter algebra (^4, R), we define on the K- vector space underlying A the following 
so-called double Rota-BcLxter product 

a *ij 6 := R{a)b + aR{h) - ah, \/ a,h e A. (3) 



Matrix Representation of Renormalization in pQFT, August 20, 2005 



5 



Then the vector space A equipped with the product *ji and operator R is again a Rota-Baxter 
algebra (of weight one), denoted by Ar and called the double Rota-Bctxter algebra of A. 
Further, The Rota-Baxter map R becomes an algebra homomorphism from Ar to A: 

R{a*Rb) = R{a)R{b), Wa,b e A. (4) 

For the Rota-Baxter map R we find R{a *r b) — —R{a)R{b). 

2.2 Complete Rota-Bcixter algebrsis and Spitzer's identity 

A complete filtered Rota— Baxter algebra [25] is defined to be a Rota-Baxter algebra {A, R) 
with a complete decreasing filtration of Rota-Baxter ideals {An}n>o- So we have 

*^n+l ^ Am A-mA-n ^ A.m-\~ni -^(*^n) ^ Am ^ilillA/An — A- 

The last equation is equivalent to saying that A is complete with respect to the topology on A 
defined by the ideals An- 

Examples 1. Let (^4, R) be a Rota-Baxter algebra with Rota-Baxter map R. We have the following 
complete Rota-Baxter algebras. 

1. The power series algebra A := ^[[a^]] where the filtration is given by the degree in x and 
the Rota-Baxter operator IZ : A ^ A acts on a power series via R through the coefficients, 

^ (En>0 ■= En>0 R{(^n)x''] 

2. Let Hj^he the connected graded Hopf algebra of Feynman graphs. The regularized functionals 
Hom(7i^_;r, ^4) where the filtration is given by the grading oiHjr and the Rota-Baxter operator 
acts on a hnear map / : Hjr — > A by acting on the target space image of /. See Theorem 16; 

3. The upper triangular matrices with entries in A where the filtration is given by the number 
of zero subdiagonals of the matrices and the Rota-Baxter operator acts on a matrix entry by 
entry. See § 2.3. 

By the completeness of {A, R) , the functions 

exp : ^1 — > 1 + exp(a) 

log: 1 + Ai^Ai, log(l + a) : 
are well-defined and are inverse of each other. 



n=0 



■~ ^ n! ' 

{-ay 



n=l 



n 
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The algebraic formulation [3, 13, 33, 45] of the classical Spitzer identity [47] was first given for 
a commutative Rota-Baxter algebra A. In terms of the complete filtered commutative Rota-Baxter 
algebra A — A[[x]\ in the first example above, it takes the form 

oo 

exp ( - 7^(log(l + ax))) = ^(-l)"(i?a)["]x". Vox e Ai. (5) 

n=0 

Here we inductively define 

(i?a)["+^l := i?((i?a)M a) 

with the convention that (i?a)[°l = 1. We note that the element b := E~=o(-l)"(-^«)'"'^" ^ 
uniquely solves the recursive equation 6 = 1 — TZ{bax). 

If {A, R) is a complete filtered commutative Rota-Baxter algebra of weight zero, i.e., R fulfills 

R{a)R{b) = R{R{a)b + aR{b)) , Va, 6 e A, (6) 

then Spitzer 's identity (5) reduces to the identity 

oo 

exp ( - R{a)x) = ^(-l)"(i?a)["]x", Va e ^, 

n=0 

which is well-known in the context of linear differential equations, where R is the Riemann integral. 
It follows naturally from (6) since R{a)" — n!(i?a)["l for a E A. 

In the following theorem, we generalize the Spitzer identity to non-commutative complete filtered 
Rota-Baxter algebras. The essential difference with the commutative case is the map x : — > 
appearing inside the exponential. This map was introduced in [24] and is defined recursively by 

Xiu) ■.= u- BCH{Rixiu)), R (x(w))) , Vw G A (7) 

using the BcLker-Campbell-HausdorfF (BCH) formula 

exp(a;) exp{y) = exp (^x + y + BCH{x, y)) 

which is a power series in x,y e Ai of degree 2. One finds a simpler recursion for x, using the 
factorization property implied by the x nieLp on A. 

Lemma 2. [24] Let A he a complete filtered K-algebra. K : A ^ A is a linear map. The map x 
in (7) solves the following recursion 

Xiu):^u + BCH{-Kixiu)),u)), u e A,. (8) 

Proof. In general for any u E Awe can write u = K{u) + (id_4 — K){u) using linearity of K. The 
map X then implies for u E Ai that exp(M) = exp [K{x{u))) exp Further, 

exp (X(x(ii))) = exp ( - K(xH)) exp(ii) 
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= exp ( - K{x{u)) +U + BCH{-K{xiu)),u)) . 
Bijectivity of exp map then implies, that 

X{u)-Kix{u)) = -Kixiu))+u + BCH{-K{x{u)),u). 
Prom which Equation (8) follows. □ 

More details on this recursion can also be found in [25, 40]. The factorization in the complete 
filtered algebra A follows from the map %. Replacing the linear map X by a Rota-Baxter operator 
R, we arrive at Spitzer's identity for non-commutative Rota— BcLxter algebra. 

Theorem 3. [25] Let {A, R, An) be a complete filtered Rota-Baxter algebra of weight one. Let 
a & Ai and recall that R :— id^ — R. 

1. The equation h — 1 — R{ha) has a unique solution b — exp ( — i?(x(log(l + a)))) . 

2. The equation b' — 1 — R{ab') has a unique solution b' — exp ( — i?(x(log(l + a)))). 

When {A, R) is commutative, the map x reduces to the identity map, giving back Spitzer's 
classical identity (5). In general we have Atkinson's theorem [2], giving a decomposition. 

Theorem 4. [2, 25] Let {A, R) be a complete filtered Rota-Baxter algebra. For solutions b and b' 
in items (1) and (2) of Theorem 3, we have 

6(1 + a)b' = 1, that is, (1 + a) = b'^b''^. (9) 

// R is idempotent: R^ = R, then this is the unique decomposition of 1 + a into a product of an 
element in 1 + R{Ai) with an element in 1 + R{Ai). 

Recall that R :— id^—R is a Rota-Baxter operator if and only if R is. Further R :— id^—R — R. 
Thus by exchanging R and R in the definition (7) of x and in Theorem 3, we have the following 
variation which will be useful later. 

Theorem 5. Let {A,R,An) be a complete filtered Rota-Baxter algebra of weight one. Define 
X '■ Ai ^ Ai by the recursion 

X{u) -.^u- BCH{R{x{u)),R{x{u))), W e A^. (10) 

Let a & Ai. 

1. The equation b — 1 — R{ab) has a unique solution b — exp ( — i?(x(log(l + a)))). 

2. The equation b' — 1 — R{b'a) has a unique solution b' — exp ( — i?(x(log(l + a)))). 
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3. For solutions b and b' in item (1) and (2), we have 

b'{l + a)b = 1, that is, (1 + a) = 6'"^6"\ (11) 

When R is idempotent, this gives the unique decomposition of 1 + a into a product of an 
element in 1 + ^(^i) with an element in 1 + R{Ai). 

The two decompositions in Theorem 4 and Theorem 5 are simply related as follows. 

Proposition 6. Let .4.°^ be the opposite algebra of A, with product defined by a . b := ba. Let 

O : ^°P,0(a) = a, be the canonical antihomomorphism of Rota- Baxter algebras. For a G Ai, 
let 1 + a = b'^^b'~^ be the decomposition in Theorem 4, with b (resp. b' ) being solution of item (1) 
(resp. (2)) of Theorem 3. Then in 

1 + a = 0(1 + a) = 0{b'-^) .0{b-') = 0{b')-' ■ 0{b)-\ 

the factor 0{b) (resp. 0(b')) is the solution, in the opposite Rota-Baxter algebra (A'^'^jR), from 
item (2) (resp. item (1)) of Theorem 5. 

Proof. We just need to note that, under the anti-isomorphism O : A^ the defining equations 
of X, b and b' in Theorem 3, with multiplication in A, are sent to the defining equations of x, b and 
b' in Theorem 5, with multiplication in □ 

We give another variation of Theorem 3. 

Proposition 7. Let b and b' be solutions of the equations in item (1) respectively (2) of Theorem 
3. Let a = (1 + a)"^ - 1. Then 

1. b~^ = exp (i?(x(log(l + a)))) is the unique solution of the equation c = 1 — i?(ac) . 

2. b'~^ — exp (i?(x(log(l + a)))) is the unique solution of the equation d — 1 — R{c'a). 

3. Further 

b-^ = 1 + R{ab'), b'-^ = 1 + R{b a). 

Proof. (1) Since a is in Ai, the equation c—1 — R{dc) has a unique solution. We just need to check 
that the solution c is the inverse of b. Since b satisfies b — 1 — R{ba), we have 

be = {I - R{ba)){l - R{dc)) 

= 1 - R{ba) - R{dc) + R{ba)R{ac) 

= 1 - R{ba) - R{ac) + R{baR{ac)) + R{R{ba)ac) - R{baac) 

= 1 - R{ba{l - R{ac)) - R{{1 - R{ba))ac) - R{baac) 

— 1 — R{bac) — R{bdc) — R{badc) 
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= 1 — R(b{a + d + aa)cj , 



since a + d + ad = a+{l + a)d = a — a = 0, we get be — 1, as needed. 

(2) The prove of the second statement is the same. 

(3) Note that 

db-^ = (b'b - l)b-^ = b' - b-^ = 6' - (1 + a)b' = -ab'. 
So by item (1), we have 

= 1 - R{db-^) = 1 + R{ab'). 

□ 

By Theorem 5.(1), the equation c = 1 — R{dc) in Proposition 7 has a unique solution 
c = exp ( - i?(x(log(l + a)))) = exp ( - i?(x(log(l + a)-^))) = exp ( - i?(x(-log(l + a)))). 
By the bijectivity of log and exp, we have 

R{x{u)) = -R{x{-u)), V«e A. 

Thus we obtain 
Corollary 8. 

^(xN + x(-«)) = 0, V«e^i. 

2.3 Decomposition of triangular matrices 

Let Ahedi, commutative K-algebra. In the following one might replace "upper" by "lower" without 
restriction. The algebra oi n x n matrices with entries in A is denoted by M.n{A). We have the 
subalgebras M.n{A) C Ain{A), 1 < n < oo, of upper triangular matrices. We also let M.'^{A) 
denote the algebra of cxo x oo upper triangular matrices. 

The subset 1 < n < oo, of upper triangular matrices with unit diagonals, i.e. those 

a G such that an — 1, i = 1, . . . ,n, form a group under matrix multiplication. The inverse 

of q; = (aij) G DJlniA) is given by the well-known, recursively defined inversion formula for upper 
triangular n x n matrices 

{a~-^)ij = -aij - ^ {a~%kOikj- (12) 

k=i+l 

Here commutativity of the algebra A is needed. 

For each n < oo, the algebra M^iA) carries a natural decreasing filtration in terms of the 
number of zero upper subdiagonals. We denote by A4'!^{A)i the upper triangular matrices with the 
main diagonal being zero, that is, the strict upper triangular matrices. Let A4'!^{A)k, k > 1, denote 
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the ideal of strictly upper triangular matrices with zero on the main diagonal and on the first k — 1 
subdiagonals. We then have the decreasing filtration 



with 



Ml{A),Ml{A)^GMl{A),^m. 

We also have Wln{A) = 1 + A4'^{A)i, here 1 denotes the n x n unit matrix. It is easy to see that 
for any n < oo, the filtration is complete, that is, A^^(A) is complete with respect to the topology 
defined by the ideals M';,{A)k, k>0. 
Thus the maps 

(13) 



exp : Ml{A), ^ MniA), exp(Z) = 
log : MniA) ^ Ml{A)^, log(a) = - ^ L^^. 



(14) 



k=l 



are well-defined and are the inverse of each other. We denote Za = log(Q;) for a G OJt„(A). 

Now let A be a Rota-Baxter algebra with Rota-Baxter operator R. We define a Rota-Baxter 
map 71 on A4^{A) by extending the Rota-Baxter map R entrywise, i.e. for the matrix a = (ajj) G 
Ml{A), define 

7e(a) = (i?(a.,)). (15) 

Theorem 9. The triple {Mn{A),7l,{Mn{A)k}k>i) forms a complete filtered Rota-Baxter algebra. 

Proof. We only need to show the Rota-Baxter relation for TZ. For distinction, we use brackets [ 
and ] for matrix delimiters. Let a — [aij] and /3 = be in M'^{A). By the entry- wise definition 
of TZ in (15) and the Rota-Baxter relation, we have 

7^(«) n{p) = [Ria,,)] [i?(A,)] 

k 

= (i?(i?(aife)/?fe,) + R{aikR{hj)) - R{aik(5kj))\ 

k 

= [j2R{R{aik)Pkj)] + [Y^RioiikRWkj))] - [5^^(aife/3fe,) 

k k k 

= [R{J2RMf^kj)] + [R{J2(^ikR{(^kj))] - 

k k k 

^TZ^'^R{aik)Pkj +TZ^'^aikR{Pkj) - T^^y^.aikPkj^ 

k k k 

= 7^(7^(K])[A,■]) +7^(K•]7^([A,•])) -7^([a,,■] [a,]) 

= n{n{a)(3) + n{an{(3)) - n{a(3) 
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as needed. 



□ 



Then we can apply Theorems 3, 4 and 5 to obtain decompositions of upper triangular matri- 
ces. For later applications to matrix Birkhoff decomposition in renormalization, we will stress the 
variation in Theorem 5. 

Corolleiry 10. Let a be in Wln{A). 

1. There is a factorization 

a = (16) 

of a into a product of an element in 1 + 7?.(A4^(A)i) and an element in 1 + TZ(^A4'^{A)i^ 
which is unique if B? — R. 

2. The factors and «_ have the explicit expression 

a+ = exp {n{x{Z^))^, a_ = exp ( - 7^(x(^a))). (17) 
Here = log(a) and x is defined in Equation (10) in analogy to %. 

3. Further, a]^ is the unique solution to the equation 

^1- n{p'{a - 1)) (18) 
and OL- is the unique solution to the equation 

P ^l-n((a-l)P). (19) 

Equations (18) and (19) are similar to the well-known recursions in renormalization theory 
where they are called Bogohubov recursions [15, 16]. We will explore their connection in the next 
section through a matrix representation of regularized Feynman characters in renormalization in 
perturbative QFT. 

Corollary 10 suggests that q;_ and q;+ in (16) can be calculated either by their exponential 
formulae (17) or directly from their recursive equations (18) and (19). We will first describe the 
recursive method to find the factor matrices q;_ and «+. Later in subsection 3.5.1 we will calculate 
the factor matrices using the BCH-recursion x in (10). 

As an example we first consider a straightforward 2x2 factorization of the matrix a e 9Jt2(^), 
which simply follows from R + R = id^. 
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The case of 3 x 3 matrices is already more telling. For a given a 
I a b\ 



/ 1 a b\ 

1 c 
\0 1/ 



in MsiA), the 



equation that «_ := | 1 c 

1/ 




should fulfill is (19): 

// a b\ 



b\ 


/ 1 


a 


c 





1 


0/ 


Vo 






Equating the two sides entry-wise, we simply obtain, 

a = -R{a), c = -R{c), b = -R{ac + 6) = R{aR{c)) - R{b). 

So 



/l -R{a) -R{b) + R{aR{c)) \ 







\ 

and by using Equation (12), we obtain 

fl R{a) R(b) + R(aR{c)) - R{a)R(c) \ 
1 R{c) 

1 



\ 
We similarly use (18) to find 



-R(c) 
1 



/ 1 R{a) R{b) + R{R{a)c) - R{ac) 

1 R(c) 
\ 1 



/ 1 -R{a) -R{b) + R{R{a)c) \ 



\0 



-R{c) 
1 



and then use (12) to find 



fl R{a) R{b) + R{aR{c)) - R{ac) 



a+= 1 R{c) 
\ 1 

We thus obtain the unique factorization in (16): 

fl R{a) R{b) + R{aR{c)) ~ R{ac) \ 
1 R{c} 

Vo 1 ; 



a = a+a 



-1 



/ 1 R{a) R{b) + R{R{a)c) - R{ac) \ 

1 R{c) 
\ 1 

(20) 



To compute a± effectively in general, we have 
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Theorem 11. Let a e QJt„(>l). 

1. The equation /? = 1— 7^((q; — l)/3) has a unique solution (3 = {J3ij) where 



3-^ 



4- = -R{a,,) -J2 E {-l)''+'R{auMo^hi,R{o^i,is ■ ■ ■ R{ai,.d) •••)))• (21) 

k=2 i<li<l2<---<lk-i<j 

2. The equation /9' = 1 - {a - 1)) has a unique solution (5' — where 

j-i 

k=2 i<h<h<-<lk-i<j 

Proof. We will prove Equation (21). The proof for the second equation is the same. Comparing 
two sides of /9 = 1 — 7?.((q; — l)/9), we have 

{0, i > j, 

1, i = j, 

-R(aij + Ei<„<:,- OiiuPuj) , i<j 

So we just need to prove Equation (21) for i < j. For this we use induction on j — i > 1. There is 
nothing to prove when j — i — 1. Assuming the equation holds for j — i < m, then for j — i — m + 1, 
we have 

i<u<j 

^-R{a,j) -R(j2^i-{- J2 (-i)''"'^K/i^(«/i/. • --Ri^^h-d) ■ ■ •))))■ 

\<u<j k=2 u<h<l2<-<lk-i<3 

This is what we need. □ 

The same formulae hold when R and R are exchanged. These formulae will be useful in the 
next section when we consider anti-homomorphisms of Rota-Baxter algebras. Let us also remark 
at this stage that the simplicity of the formulae in Theorem 11 suggests a computational advantage 
of working with matrices over working directly with regularized Feynman characters. 



3 Matrix calculus in perturbative QFT 

The perturbative approach to quantum field theory provides theoretical predictions which match the 
experimental data with very high precision. At the heart of this scheme lies the idea to approximate 
physical quantities by power series expansions in terms of a supposed to be small parameter called 
the coupling constant, which measures the strength of interaction in the physical system under 
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investigation. Of course, this idea raises immediately the question of convergence of such power 
series, but for such issues we refer the interested reader to more sophisticated presentations of 
perturbative QFT. The power series expansion is organized by Feynman graphs which serve to 
labels the terms in the power series. A Feynman graph is a collection of vertices and edges, reflecting 
the interaction and propagation of particles, respectively. For each term in the power series the 
number of loops of the corresponding Feynman graph fixes the order in the power series. The terms 
of the power series, called Feynman amplitudes, follow from Feynman rules which provide a way to 
translate a Feynman graph into a Feynman integral. 

In general this perturbative ansatz is plagued with so-called ultraviolet divergencies, i.e. the 
Feynman integrals corresponding to Feynman graphs diverge in the limit of large momenta or 
equivalently small distances, and therefore seem to be useless in physics. For now we ignore the 
infrared problems appearing in theories with masslcss particles. 

Renormalization theory in perturbative QFT was developed to cure these divergencies in a 
meaningful way. Actually, it consists of two steps, first a regularization prescription meant to 
control the divergences by extending the target space of the Feynman rules from the base field C 
to a particular algebra of regularized Feynman amplitudes. As a main example we mention here 
dimensional regularization, where this algebra is that of the Laurent series. Second, on the new 
target space algebra, a renormalization scheme is introduced, that isolates the problematic pieces, 
i.e., the divergent part. The renormalization process itself is made up in a recursive manner based 
on the self similar structure of Feynman graphs. Hereby we mean the fact that graphs of lower 
order -in terms of the number of loops- appear inside Feynman graphs of higher order. 

Unfortunately, despite its impressive successes, renormalization was stigmatized, especially for 
its lack of a sound mathematical underpinning. One reason for this weakness might have been the 
fact that Feynman graphs in itself appeared to be unrelated to any mathematical structure that 
may possibly underlie the renormalization prescription. 

This changed to a great extend with the original paper by Kreimer [35] , followed by the work of 
Connes and Kreimer [16, 17, 18, 19]. The combinatorial- algebraic side of the process for renorma- 
lization is captured via a combinatorial Hopf algebra of Feynman graphs, essentially characterized 
by the coproduct map which organizes the decomposition of a Feynman graph into its subgraphs in 
a firm way. The analytical side of renormalization, i.e. the Feynman rules providing the Feynman 
integrals, is imbedded as algebra homomorphisms in the dual space of this Feynman graph Hopf 
algebra. 

3.1 Connes— Kreimer renormalization Hopf algebra 

Let us briefly summarize the main results of the Hopf algebra description of renormalization theory. 
For a general introduction to Hopf algebras the reader might want to consult the standard references 
like for instance [1, 29, 48]. Details about Hopf algebras in renormalization theory can be found in 
the survey articles [26, 40]. 
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3.1.1 Feynman graphs and decorated rooted trees 

Let us mention for completeness how Feynman graphs and decorated rooted trees are related. We 
assume that the reader is somewhat familiar with the former. The latter will be introduced in the 
next section. For more details on Feynman graphs, such as the notion of ultraviolet (UV) subgraphs, 
we refer to the standard literature such as [15]. Kreimer [37], and then Connes-Kreimer [16, 17] 
well-developed this link. 

A Feynman graph is a collection of internal and external lines, or edges, and vertices of several 
types. A proper subgraph of a graph is determined by proper subsets of the set of internal edges 
and vertices. Of vital importance are so-called one particle irreducible (IPI) Feynman graphs, a 
connected graph that cannot be made disconnected by removing one of its internal edges. In general, 
a Feynman graph F beyond one-loop order, is characterized essentially by the appearance of its UV 
Feynman subgraphs ji C F. For instance, two proper Feynman subgraphs 71, 72 C F might be 
nested, 71 C 72 C F, or disjoint, 71 fl 72 = 0. This hierarchy in which subgraphs appear inside 
another subgraph is best represented by a decorated rooted tree. The example below is borrowed 
from (^^^.^-theory. 




Figure 1: Example of a 3- loop graph with two nested disjoint UV subgraphs from lyO^^ -^-theory. 

The two UV subgraphs Q (little boxes inside) are disjoint, but nested inside another Feynman 
graph of the same type, >0<- The rooted tree on the right represents this hierarchy. There is a third 
possibility, that subgraphs might be overlapping. Such Feynman graphs are represented by linear 
combinations of decorated rooted trees. Let us illustrate this with an example from yjg^-^-theory. 




Figure 2: 2-loop example with overlapping UV subgraphs from ^jg^^^-theory. 

For a detailed treatment of this special and important case of overlapping structures we refer to [37]. 
The essential combinatorial operation on the set JF of (equivalence classes of) Feynman graphs in 
the process of renormalization is a particular decomposition of such a Feynman graph into its UV 
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subgraphs, well-known to the practitioners under the name of Bogoliubov's R-operation, or its 
solution by Zimmermann's forrest formula [50]. The concept of the combinatorial Hopf algebra of 
Feynman graphs, TYjf, enters via its coproduct map, denoted A : Tijr Tijr ® Tijc, which organizes 
such a decomposition in a mathematical sound way, e.g. for the last examples from y^g^-^-theory 
we find 



A 



which must be compared with Bogoliubov's formula for the counter term C(r) of a Feynman graph 
r. Only here we apply this notation for the counter term used in [16], and we use a symbolic graph 
notation. R denotes the renormalization scheme map 





+ R\R 




+ R 




Kreimer [35] was the first to realize in renormalization this underlying Hopf algebra structure. 

From a conceptual point of view, and having the non-specialist in mind, we feel that it is useful 
to work with decorated rooted trees. Also, this underlines the generality of our results, applicable to 
any theory in perturbative QFT and its renormalization by choosing the particular set of decorations 
dictated by the theory itself. Nevertheless, in terms of applications of Connes-Kreimer's Hopf 
algebra of renormalization to physics it is most naturally to formulate the Hopf algebra directly 
on Feynman graphs, once a theory has been specified. We included a simple calculation in section 
3.4.3 using an example from y^^^-^-theory. Also, we refer the reader to the companion paper [23] 
with J. M. Gracia-Bondia and J. C. Varilly, which contains detailed applications of the presented 
results to renormalization in perturbative QFT. We should warn the reader, that in [23] a slightly 
different notation is used. 



3.1.2 Connes— Kreimer Hopf algebra of rooted trees 

A rooted tree t is made out of vertices and nonintersecting oriented edge, such that all but one 
vertex have exactly one incoming edge. We denote the set of vertices and edges of a rooted tree 
by V{t), E{t) respectively. The root is the only vertex with no incoming line. The empty tree is 
denoted by Ir- Each rooted tree is a representant of an isomorphism class, and the K- vector space 
freely generated by the set of all isomorphism classes will be denoted by T. 



Ir . I 



i A 1 1^ A A n 1 ... •^ ... 



Definition 1. The K-algebra of non-planar rooted trees, denoted by Tir, is the polynomial algebra, 
generated by the symbols t, each representing an isomorphism class in T. The unit is the empty 
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tree Ir- The disjoint union of rooted trees serves as a product, denoted by juxtaposition, i.e. 
mArit',t) =: ft. 



There exists a natural grading on Tij- in terms of the number of vertices of a rooted tree, 
:= On forests of rooted trees, we extend it to #(ti---tn) '■= Yl^=iH^i'^i)- ^^^^ 

Tir = IK © ©n>o becomes a graded commutative K-algebra. 

To define a coproduct on Tir-, we introduce the notion of admissible cuts on a rooted tree. A cut 
Ct on a rooted tree t G T is a subset of edges q C E{t). It becomes an admissible cut if and only 
if along a path from the root to any of the leaves of the tree t, one encounters at most one element 
of Q. By removing the subset q of edges, each admissible cut q produces a forrest of pruned trees, 
denoted by The remaining part, which is a single rooted tree linked to the original root, is 



denoted by Rc^.. We exclude the cases, where q 



), such that i?ct 



and the full cut, 



such that Rr 



Pc^ = t. The rooted tree algebra Ht is equipped with a bialgebra structure by 



defining the coproduct A : Ht 'Hr ® 'Hr in terms of all admissible cuts Ct of a rooted tree t: 



A(t) = t ® Ir + Ir ® t + ^Pc, 



R. 



Cf 



(23) 



ct€Ct 



For example. 



A(.) 

A(A) 
A 



.O Ir + Ir ® . 

I O Ir + Ir ® I + . O . 

[(8)lr + lr®S + .®I + I®. 
A® Ir + Ir® A + 2.®I + 

+ .® j + lol 



(24) 



Ir + 1 



r 



Latter we will use the shortened notation A(t) = J2{t) ^(i) ® ^(2)- For products of trees we demand 
the compatibility A(ti . . . t„) := A(ti) . . . A(t„). 

Remark 12. It is important to notice that the right hand side of A(t) G Tir ® ?ir linear. 
Therefore we can write, for the K-vector space KT spanned by T, 



The counit e : Hq — > K simply maps the empty tree Ir to Ir G IK and the rest to zero 



(25) 




(26) 
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Theorem 13. [16, 35] The algebra Ht equipped with the above defined compatible coproduct A : 
Hq — > Ti.r^'Hr and counit e : Hq — > IK forms a connected, graded, commutative, non-cocommutative 
bialgebra. In fact it is a Hopf algebra, with antipode S : Ht — 'Hr defined recursively by S{lq-) — It; 

s(t)^-t-j:^^^c,s(Pc,)Rc,. 

We denote the restricted dual of Hr by H^. It contains linear maps from Hr into K. Let A be 
a commutative K-algebra, we also consider the linear functional Hom(?ir, A). Equipped with the 
convolution product: 

/ * ^ := ruAif ®g)/^:Hr ^Hr^Hr ^ A®A^ A (27) 

Hom(7ir, ^) becomes a non-commutative EC-algebra with unit given by the counit e of Tir- 

A hnear map : Hr — > ^ is called a chciracter if is an algebra homomorphism. We denote 
the set of characters by Ga '■— charATi-r- 

Proposition 14. The set of characters Ga forms a group with respect to the convolution product 
(27). The inverse of (p E charA'Hr is given by := (p o S. 

A linear map Z : Hq — > A is called a derivation, or infinitesimal chciracter if 

Z(t,t2) = ZitMh) + e(£i)Z(t2), ti eHr: i = 1, 2. (28) 

We have Z{lq-) = 0. The set of infinitesimal characters is denoted by gA '■= dchar aHt ■ For any 
Z G dchar a'Ht and t G Tir of degree #(t) = n < oo, we have Z*"^{t) = for m > n. This implies 
that the exponential cxp*(Z)(t) := Xlfc>oTr(^)' ^ ^ dchar A'Hr, is finite, ending at k = i^{t). 
Given the explicit base of rooted trees generating Tir, the A-module of derivations dchar A'Hr is 
generated by the dually defined infinitesimal characters, indexed by rooted trees 

Zt{t'):^5t,t>lK, t,t' eT. (29) 

Proposition 15. The set gA = dchar A'Hr defines a Lie algebra when equipped with the commuta- 
tor: 

[Zf, Zf] :— Zfi ic Zt'i — Ziii -k Zii 

= Y,{n{t!,f-t)-n{f,^-t))Z,, (30) 
ter 

where the n{t', t"; t) & N denote so-called section coefficients, which count the number of single 
admissible cuts, \ct\ — 1, such that Pc^ — t' and Rc^ — t" . 

Ultraviolet divergencies demand for a regularization of the theory by replacing the base field K 
in Hr as target space, by a commutative, unital Rota-Baxter K-algebra (A, R). As an example we 
mentioned earher the commutative Rota-Baxter algebra of Laurent series A :— C[e~^,£]]. 

The following theorem recalls the Birkhoff decomposition of Hopf algebra characters from the 
work of Connes and Kreimer in the context of complete filtered Rota-Baxter algebras. 
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Theorem 16. [24] Let Ti. be the connected graded Hopf algebra of rooted trees, Ht or Feynman 
graphs, Hj^, and let {A, R) be a commutative Rota-Baxter algebra with Rota-Baxter operator R 
of weight 1. Let A :— Hom{H, A). Define a decreasing filtration {An}n on A in duality to the 
increasing filtration on H by grading. Define TZ : A^ A by TZ{f) — Ro f , f E A. For a regularized 
character e Ga, let 

= 0:^*0+ (31) 

be the algebraic Birkhoff decomposition of Connes-Kreimer [17] that gives the counter term 0_ and 
renormalization . 

1. The triple {A,TZ, {Ak}k) is a complete Rota-Baxter algebra under convolution (27). 

2. The decomposition (31) is the factorization of (j) in Theorem 4, with 0_ e e + TZ{Ai) and 
(f)'^ e e-\-TZ{Ai). For idempotent Rota-Baxter map R we have uniqueness of the factorization. 

3. Explicitly we have 

(/._ = exp* ( - n{x{Wm)), r+ = exp* ( - n{x{\og*m)) 

where :— log*(0) G gA with log* and exp* defined by convolution product (27), and x is 
the BCH-recursion (7). 

4. Further, 0_ and (f)'^ solve Bogoliubov's recursions, i.e. item (1) respectively item (2) of 
Theorem 5 on Spitzer's identity, for (p — e & Ai: 

(/._ = e - 7e(0_ ★ (0 - e)) = e-:^((0-e)*0;^) 

= e + 7^(exp*-(-x(log*(0)))-e) = e - :^(exp*- (x(log*(0))) - e) 

Proposition 17. [24] Let (p e Ga be a Feynman rules character, with :— log*(0) e gA ■ The 
map 

b[ct>] := exp*^ ( - x{Z^)) - e (32) 

defined in terms of the exponential with respect to the double Rota-Baxter product represents Bo- 
goliubov's preparation map, also known under name R-operation. 

Here exp*'^ denotes the exponential map with respect to the Rota-Baxter double product (3) 
defined in the Rota-Baxter algebra (.4,7^, {Ak}k)- 

3.2 Matrix representation of linear functionals 

Recall that a Hopf algebra H is called filtered if there are K-subvector spaces n > Oof 7i such 
that 

1. 7^^") C -^("+1); 
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2 



3 



4 



p+q=n 



7i is called connected, if in addition 7i*^°^ = K. Then for any x E 7i'-"'\ we have 



A{x) := A(a;) - X(g)l-l(g)XG ©p+q=n,p>0,q>0' 



(33) 



Definition 2. A subset X of 7i is called a (left) subcoset if X is K-linearly independent and if 
KX is a left subcomodule of Ti. A subcoset X is called filtration ordered if X is given an order 
that is compatible with the order from the filtration of Ti. In other words, if i < j and Xj is in Tin, 
then Xj is in Tin- A subcoset X is called a 1-subcoset if ly^ is in X. 

It is clear that X is a left subcoset of Ti means that the coproduct A on restricts to Ax : 
X — > ® X. In other words, for any x G X, with Sweedler's notation Ax{x) = ^(^.) ^^(i) ® X(2), 
we have X(2) G KX. 

This definition reflects the combinatorial character of certain filtered Hopf algebras. Objects, 
mostly of graphical type, with specific substructures are disentangled. For an example see the 
coproduct example of the Feynman graph, ^CD^, in subsection 3.1.1. The two subgraphs ^ and 

are replaced by a vertex, ^ , with the cograph -O as a result. 

Our matrix representation applies to any connected filtered Hopf algebra H with a comodule. 
It naturally generalizes the classical construction of modules from comodules [40, 48]. 

Examples 18. Examples of such Hopf algebras with a left 1-subcoset X include 

1. the Hopf algebra Hr of rooted trees with X = T being the set of rooted trees. See Remark 12; 

2. the Hopf (sub) algebra of ladder trees TCm C TCt with X = id C T being the set of ladder 
trees, i.e. trees whose vertices (except the root vertex) have only one incoming and at most 
one outgoing edge; 

3. the Hopf algebra of Feynman graphs with X = JF being the set of one particle irreducible 
Feynman graphs. This follows from a remark similar to Remark 12; 

4. one of the above Hopf algebras with X being the first n {n > 1) elements in the subset there, 
with an ordering compatible with the filtration of the Hopf algebra. 
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3.2.1 The representation anti-homomorphism map 

Let A he a commutative K-algebra. Denote AX = A (E) KX which is a free A-module with 
basis X. For a Hopf algebra Ti with a left 1-subcoset X, we will define a linear map \1'x,a '■ 
Hom(7i, A) — > End{AX), eventually giving rise to a natural representation of Hom(7i, A) in terms 
of upper triangular matrices with entries in A. In the following the subscript X will be suppressed 
when possible. 

Definition 3. Let 7i be a Hopf algebra with a filtration ordered left 1-subcoset X. A is a com- 
mutative K-algebra. The linear map '^a,x '■ HomiTi., A) — > End{AX) is defined by taking the 
composition 



-.A^KX A^n^KX ""^^f^""^^, A0A0KX A^KX (34) 

for / e Hom(H,A). 

In particular, for x e KX with A{x) = X^(^) a;(2), we have 

*a[/](x) = (35) 

This justifies the short-hand notation f -kid^x for the composition (34) defining 

As an example, consider the Hopf algebra Hr of rooted trees with X — T. Using the notation 
in (23), we have, for / e Hom(7i-r, ^) and t eT, 

^A[/](t) = f{lr)t + f{t)lr + fiPc^)Rc, (36) 

ct&Ct 

which is in AT. 



3.2.2 Coproduct matrix for filtered Hopf algebras 

With the natural basis X oi AX, the map gives a matrix representation of Hom(7i, A). Ex- 
plicitly, fix a linear order {x,i}i>i of the left 1-subcoset X that is compatible with the filtration of 
the Hopf algebra Ti. Then the coproduct A of 7i restricted to X writes 

oo 

A{xj) = Y^Xij(^Xi (37) 
1=1 

for uniquely determined Xij in V.. Note the order of i and j. For all i < j, ^{Xy) < #(xj). This 
leads to the 

Definition 4. Let be a filtered Hopf algebra. Fix a linear order {xi}i>i of the filtration ordered 
left 1-subcoset X. We define the |X| x |X| matrix 

Mn := Mn,x := {Xij) 
in M-ixiCH), called the coproduct matrix of H (with respect to X). 



Matrix Representation of Renormalization in pQFT, August 20, 2005 



22 



Then Myi is upper triangular since Ti is filtered. For a different choice of basis, we get conjugate 
matrices. Under a fixed ordering (xi, ■ ■ ■ , a;„, ■ ■ ■ ) of X, we obtain an isomorphism 



AX ^ A 



sending x„ to the unit column vector with 1 in the n-th entry and zero elsewhere. Here we use 
the familiar bra-ket notation to denote abstract vectors as column vectors. We likewise obtain the 
isomorphism 

End(AX) ^M^x\i.^) 
sending ^^aI/] to the upper triangular matrix 



We further have an isomorphism 



Hom(AX,A) ^ A\^\ 



(38) 



sending the dual basis x*, defined by x* (x^) = 5n^mi to the row vector (x„| with 1 in the n-th entry 
and zero elsewhere. The image of / G Hom(y4X, A) under this isomorphism is denoted by (/|. So 
we have 

(/| = (/(xi),---,/(x„),---). 

We will often use these isomorphisms as identifications when there is no danger of confusion. In 
particular, / is often identified with 

As an illustrating example we consider again the Hopf algebra Tir to rooted trees with a trun- 
cated X being 



(6) 



|ei := Ir, 62 := ., 63 := I, 64 := |, 



65 := A , 66 := 



(39) 



with the displayed linear order. Then we have the corresponding column vectors 



{|lr), |.), |I>, \\), |A>, |{^} 



(40) 



In particular. 



|lr) 



/ 1 \ 





and 



(lr| = (l,0,---,0). 
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M; 



H 



( Ir 





V 



62 

Ir 







64 

es 

62 

Ir 





65 
6262 
262 



Ir 




66 
64 
63 
62 



Ir 



Prom the coproducts in (24), we obtain the truncated unipotent coproduct matrix from Defini- 
tion 4 

63 
62 
Ir 






(41) 



For an / G Hom(7i7-; ^) we then obtain its representation /(g) in terms of an upper triangular 
matrix by applying / to Myi entry by entry. In each column we just see the image of the left hand 
side of A(6j) under /. 

For an infinitesimal character Z G (7 we find particularly simple nilpotent matrices 



^\Z\{t) = Z{t)\ 



r 



ct&Ct 
\ct\=\ 



Z{t')t" 



(42) 



Here, due to relation (28) we only need to consider single admissible cuts, \ct\ = 1. The sum on 
the right hand side of (42) therefore goes over all decompositions of the tree t, resulting from the 
removal of exactly one edge. The tree t' denotes the pruned subtree of tree t, and t" the denotes 
the tree, which is still connected to the root of t. For the generators Zt we find 

^Zt,]{t2) = 6t,,t,lT + Yl ^^*i'*2^2 = 5ti,t2lr + ^n(ti,t;t2)t. 

|ct2l=l 



As an example, let us calculate 



nz.]{A) 



2Z.(.)I = 2I 
Zj(64). = . , 



^Z.]{e,)=Z.i,): 

nZ^]{e,)=Zf).= 



^Z.]{ee)=Z. 



.)Li 



^[Z^]{e,)=Z0 



I. 



For the generators Zt of the Lie algebra of derivations we find in this particular example the following 
matrices. 



( ° 














1 \ 




/o 











1 


o\ 




/ 








1 





o\ 

























































1 
























































































































































































































lo 














0) 




Vo 














0) 




\o 














0/ 



(43) 
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/ 1 \ 

1 

1 





Voooooo/ 



/ 1 \ 

1 

1 2 

1 



voooooo/ 



(44) 



For a character G Hom(7ir, A), i.e. an algebra homomorphism from Tir to A, 0(lr) 
Ik, we find the upper triangular matrix with unit diagonal G WIq{A) C A4q{A) 



( Ik 





V 



le2j 
Ik 







^(es) 0(64) 0(65) ^(ee) \ 

0(62) 0(63) 0(62)^ 0(64) 

Ik 0(62) 20(62) 0(63) 

Ik 0(62) 
Ik 
Ik 



(45) 



The counit, which is a character by definition (26), is represented by the unit matrix 'e 
structure of the representation matrices shows that the row vector 

(lr|0= (1,0, 0,0, 0,0)0 



1. The 



IS 



(01 := (Ik, 0(62),- ■■ ,0(e6)). 

Returning to the general case, where 7i contains a (filtration ordered left) 1-subcoset X, and 
(A, K) is a Rota-Baxter algebra, we have the following key property of ^ a,x that establishes 
the connection between the Rota-Baxter algebras Hom(7i, A) and hA^i^A). Remember that we 
suppressed the subindex X if there is no danger of confusion. 

Theorem 19. Let Ti. be a connected filtered Hopf algebra with a filtration ordered (left,) 1-subcoset 
X dTi. Let A be a Rota-Baxter algebra (of weight 1). 

1. The linear map 

■■= : Hom(7^, A) ^ M\x\{A) 

is an anti-homomorphism of Rota-Baxter algebras that is continuous with respect to the topolo- 
gies defined by the filtrations on the filtered Rota-Baxter algebras. More precisely, for any 
m > 1, there is N > 1, such that for all k > N and f G Hom(?i, A) with fiTV') = 0, we have 

^A[f]eM];^^{A)m. 



2. The first row off := ^ A,x[f] ^ M^x\iA) is (/| := (Ul / = {f{x)).^x. 
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3. If X is a generating set of the algebra 7i, then the map restricts to an injective map 

from the multiplicative group Ga — char^Ti^ of algebra homomorphism H. ^ A to End(AX). 

Proof. (1) We denote the composition in End(74X) by concatenation. Let f,g & Hom(7i^, ^4), 
and X & X. We will use Sweedler's notation of ^{x) — X{i) ® X{2) and the short hand notation 
= /*idAx in (35). By the definition of ^ a, the coassociativity of the coproduct A and the 
definition of the convolution product we have 

^> A[f]^ A[g]{x) = {f ^idAx){{g^idAx){x)) 

= {f -kidAx{9{X(^i))X(2)) 

= g{x^l)){f{Xi2)(l))X(2)(2)) 

= (5(2^(1)(1))/(2;(1)(2)))2;(2) (46) 
= f){.X{^i))Xi2) 

= ((^*/)*idAx)(2;) 

= -^Aig^fMx) 

This proves the anti-homomorphism property. 

Further, TZ is the Rota-Baxter map on Hom(?i, A), and we let TZ denote the Rota-Baxter 
operator on A1|^|(>1) by acting entry- wise. Let x ^ X with A(a;) = ® X{p,y We have, for 
/ e Hom(H,^), 

^A,x[7^(/)](a;)= W)(x(i))x(2) 
= i?(/(a;(i)))x(2) 
= :^(/(a;(i))x(2)) 

= ^(*A,X[/])(X). 

This proves that a,x is compatible with the Rota-Baxter operators. 

The continuity is verified using the fact that the linear order on X is compatible with the 
filtration of T-L. 

(2) Since X is assumed to contain 1 — Ix and any linear order of X is assumed to be compatible 
with the filtration oiH, we have 

X = {xi := 1,^2, •• •) 
for any choice of such linear ordering. So by Equation (33), we have 

A(xj) — Xj <S>1 + Xij (8) Xi- 

i>2 

Therefore the first row vector of the coproduct matrix My, is (xi, X2, • • • )■ So the first row of / is 

{fix,),f(x2),---). 
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(3) Suppose = 0. Then / = 0. By item (2), f{x) = for all xeX. Since f : H ^ A is 

an algebra homomorphism and X is an algebra generating set oi H, we have f — 0- This proves 
the injectivity. □ 



3.3 Renormalization of Feynman rules matrices 

We now consider the Hopf algebra H^r of Feynman graphs. Let be the set of (equivalence classes 
of) Feynman graphs with a fixed hnear ordering Fi := 1, • • • compatible with the filtration H^^^ 
of TijF- Then with X being either or the subset J^n of the first n-elements in the results of 
last subsection apply. 

With the representation of a regularized Feynman character : Tijr — > by an upper triangular 
matrix obtained in Theorem 19, we can apply the complete filtered Rota-Baxter algebra structure 
on the matrices to decompose the representation 0, via Corollary 10, into the inverse of the counter 
term matrix 0_ and the renormalized matrix 0+, giving rise to an analog of Connes-Kreimer's 
algebraic Birkhoff decomposition for regularized Feynman characters. Then by Theorem 19.3, these 
later matrices recover the counter term and renormalization of the regularized Feynman character 0. 

In hght of the anti-homomorphism in Theorem 19, we expect that exchanges the order 
of the Birkhoff decomposition of in (31). The precise relation between these two decompositions 
is provided by the following theorem by applying Theorem 19 and Proposition 6. We will use TZ to 
denote both, the Rota-Baxter map of Hom(7i^, ^4) and A4'!^{A), since from the context it will be 
clear which one is used. 

Theorem 20. Let (f) be a regularized character from Tij: to {A,R) and = log*(0) e gA- For a 
fixed 1 < n < oo, let X — J^n be the first n Feynman graphs and let e 3Jln(A) C A15^(A) be the 
upper triangular matrix representation of (f> given by (38). Let 



he the Birkhoff decomposition of (p E Ga in (31). Let (0_) and (0+) he the matrix representations 
of (p- G Ga respectively 0+ G Ga- 

1. Then (0_) is the unique solution of the equation 



(/)+ 



/3 = l-7^((0 



(47) 



explicitly given by 



(0_) = exp (-7^(x(^<^)) , 



and (0+) 



is the unique solution of the equation 



/3' = l-7^(/5'(0-l)) 



(48) 
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explicitly given by 

(^)~' = exp (^-n{x{z^)') . 

Here x is the BCH-recursion defined in equation (10). 

2. (f) factorizes as 

(49) 

3. The first row vector of (0-i-) is (0±(r))pg^ . This can be summarized by the linear renorma- 
lization matrix-vector equation for := (0+) and {<j)\ := 

(0+1 = (01 (^T), (50) 

following from the matrix Birkhoff decomposition in (49). 

Remcirk 21. The renormalized matrix (0+) uniquely solves the equation 

c = i-n{i^-'-i)c), (51) 

where we used Propositions 6 and 7. Instead of inverting (0+) — > (0_|_) at the end of calculation, 
here we first invert the Feynman rules character matrix 0, which we can do simply and efficiently 
using the inverse coproduct matrix, see below. 

3.4 Summciry/ Algorithm and examples 

We now discuss applications of our matrix approach to renormalization in Theorem 20. One ad- 
vantage of the matrix approach lies in the efficiency to calculate the matrices (0_) and {4>+) of the 
matrix Birkhoff factorization directly in the matrix algebra, see Corollary 10, using the recursions 
of Theorem 11. 

We outline an algorithm which is followed by examples. 

3.4.1 Calculation of the coproduct matrix and its inverse 

Let X := T he the set of rooted trees (or Feynman graphs) and fix an ordering of T that is 
compatible with the grading of T, ti := Ir, ^2, • " " ,tn,' ' ' ■ Then the coproduct A on Ht is given by 

A(i) = ^t'0i", 
it) 

or in our ordering, 

j 

1=1 
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where T^j is in Hr-, U is in T, see Remark 12, and #{Tij) < Define 



Mh = {Tij) 



to be the coproduct matrix of Hr- M-^ is an oo x oo upper triangular matrix with entries in Hr 
and unit on the diagonal. The entries Tij can be obtained by the S+ operator of Connes and 
Kreimer providing a recursive way to calculate the coproduct. One can also consider the truncated 
coproduct matrix by restriction to the subspace of the first n graphs ti, - ■ ■ ,tn and obtain a finite 
upper triangular matrix in A4'^{TC). 

The direct calculation of the renormalization matrix in Theorem 20 respectively Remark 21 
demands an inversion of character matrices. In Proposition 14 we stated the fact that the inverse 
of a Hopf algebra character e G is given by composition with the antipode anti-homomorphism, 
0-1 — (f)o S. The representation matrix of loosely speaking follows from (38), cf) :— (f)o M^. As 
we want to calculate the inverse coproduct matrix, M^"*^ = (^7^)' ^^^^ ^^^^ :— (f) o M^^, we 
may take another look at the coproduct matrix M^. 

For this purpose we go back to Definition 3, where we define the commutative K-algebra A-.— li,. 
For a filtered Hopf algebra H with a filtration ordered (left) 1-subcoset X this provides us with an 
upper triangular matrix representation of Hom(7Y,7i^). The coproduct matrix thereby represents 
the identity homomorphism id-^ : — > 



We suppressed the 1-subcoset X C Ti in the notation. 

As an upper triangular matrix, the inverse of M-j-i follows immediately in a recursive manner 
from Equation (12). But we want to make a little detour, using the aforementioned representation 
point of view in (53). This will provide us with a non- recursive simple formula for calculating M^. 

The antipode 5" in a Hopf algebra Ti is an anti-homomorphism, defined as the solution of the 
equation S * idn — e — id-^ ★ S. From a convolution product point of view the antipode is the 
inverse of the identity map, S — id^^. For connected filtered Hopf algebras the identity can simply 
be written as id-^ = exp*(log*(id7^)), using the bijectivity of exp* and log*, and therefore we have 
S — exp*(— log*(id-^)). We already mentioned the trivial fact that every algebra is a Rota-Baxter 
algebra, with Rota-Baxter operator pair id and id = 0. The space End(Ti) is equipped with two 
products, composition and convolution, both forming an associative algebra. Spitzer's identity for 
the non-commutative algebra (End(Ti),T*r) with Rota-Baxter map id-^ : End(Ti) End(Ti) then 
implies for the antipode 



(53) 



S = exp* (-id^(log*(e + (id^ - e)))) 
to be a solution of the equation b — e — id-^(6* (id-^ ~ ^))- Explicitly 



5" = e - (id^ - e) (id^^^ - e) ★ (id^ - e) - (id^^ - e) ★ (id^ - e) ★ (id^ - e) H 



(54) 
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The BCH-recursion x (7) obviously does not enter in this particular case as id-^ = 0. See [27] for 
the geometric series ansatz. Applying the matrix representation anti-homomorphism ^f-^, we find 
the inverse coproduct matrix 

^n[S] = i^ = M^i (55) 
= l-{Mn-l) + {Mn-l){Mn-l)-{Mn-l){Mn-l){Mn-l) + --- (56) 
= l-5^(M„-l)^ (57) 

n>0 

The inverse character matrix of e G therefore is given by 

0-^ = o M-^ = 1 - J](0 - 1)-. (58) 

m>0 

For a truncated coproduct matrix e M.'^{T-L), n < oo, the series on the right breaks up at order 
n, {(f) — 1)" = 0, as — 1 e A4n{A)i is nilpotent. In components, the formula for (f)~^ is 

k=l i<h<l2<-<lk<3 

3.4.2 Matrix renormalization by factorization 

Now let : Tiq > A denote a regularized Feynman rules character with image in a Rota-Baxter 

algebra [A, R). Applying to gives the Feynman rules matrix 

0:= 0(MH) = (0(r,,-))- 
Let j3 be the unique solution of the recursion 

/3 = l-7^((0-l)/3), 

as in Theorem 20. The matrix j3 can be effectively computed by Theorem 11.1. The first row vector 
of (3 is the counter term vector (0„(ti), ■ ■ ■ , (f)^{tn), ■ ■ ■)■ Then the first row vector of the matrix 
product /3 is the renormalization vector, i.e., we have the linear renormalization matrix- vector 
equation (50) of item (3) in Theorem (20) 

(lr|0/3 = (0| /3 = (0+(ti),-- - ,0+(y,---)- 
Alternatively, let be the unique solution of the recursion 

(3' ^l-n{(3' 

as in Theorem 20, again effectively computable by Theorem 11.2. Then find P'"^ which can be 
computed by Equation (58) (or recursively (12)). Then the first row of is the renormalization 
vector 

(lr|/3'-' = (0H-| = (0+(ti),---,0+(*n),---)- 
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Equivalently, using the inverse coproduct matrix to calculate the inverse Feynman rules character 

:=0(M^i) = (0(7;t1)), 

we find directly the renormalized character matrix as solution of the equation c = 1— "^((0^^ — l) c) 
of Remark 21. 

3.4.3 Examples in (^|^^^-theory 

This subsection serves to show how the above Hopf algebra consideration nicely applies to standard 
Feynman graph calculations of perturbative renormalization. In [17] Connes and Kreimer showed 
in full generality that the set of Feynman graphs T for any perturbatively treated QFT can be 
made into a combinatorial Hopf algebra Hjc of the above type. 

We will use a simplified version of (^^-theory in four dimensions as our Feynman graph toy model 
physics theory. A more detailed and refined treatment can be found in the companion paper [23]. 

The Feynman graph Hopf algebra is denoted by TYjf. As regular ization scheme we choose 
dimensional regularization. So that the space of linear functionals Hom(7ijr, A) contains maps into 
the commutative Rota-Baxter algebra A := C[£:^^,£:]] with Rota-Baxter map Rms- 

We work up to 3-loop order, by taking X to be the set of graphs 

-^(4) := {ei := |M , ■= lO) , 63 := | Y >, e, := l^^)} (60) 



identified with the corresponding column vectors. The graph is obtained by substituting a 
wine-cup Y i^to the divergent 1-loop graph xQx- 

The Feynman graphs we consider here are of purely combinatorial type in the sense that the 
external legs are not decorated by external structure, such as external momenta, spin indices, etc. 
This frees us from symmetry consideration, which otherwise would demand a bigger representation 
space JF. 

The Feynman diagram >0< has a primitive divergence and would correspond to the one vertex 
tree • decorated by this graph, the wine-cup diagram Y contains exactly one nested subdivergence 
of type Ox and corresponds to the ladder graph of length 2, I, where the root and leaf both are 
decorated by the graph xQx. The coproducts of these two graphs therefore are given in analogous 
forms to the first two expressions in (24). The graph ''^X contains three nested subdivergences and 
correspondence to the ladder graph of length 3, with each vertex decorated by O, its coproduct 
is given by 

Let G Ga denote the Feynman rules for four dimensional ip'^-theorj in dimensional regularization, 
together with minimal subtraction scheme, i.e. Rota-Baxter map R := Rms- The corresponding 
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coproduct matrix Mj^^^^, respectively character matrix = are given by 



o M- 



(4) 



/I o y ^\ 

1 o< Y 

1 o 

V 1 / 



{ 1 0(0) 0(y) ^ 
1 (/)(o) 0(y) 

1 0(0) 

\ 1 / 



This matrix has the special property that along any subdiagonal we find the same entry. This 
situation always appears if (and only if) we deal with strictly nested diagrams only. In the picture 
of decorated rooted trees this corresponds to the Hopf subalgebra of ladder trees. 
Using the counter term recursion (47) of Theorem 20 for 0, 

/3=l-7e((0-l)/3) 

and applying formula (21) for its solution, we find 



a-.-- 



where 



/ 1 


-i?((/)(0<)) 


/5l3 /3l4 \ 





1 


-i?(0(O)) /?24 








1 -i?((/)(0<)) 


Vo 





1 / 


/5l3 


=-^(0(y: 


) +i?(^0(>O)i?(0(o))), 






/324 = /3l3, 



and 



/3i4 = -i?(^0('^)^+i?(^0(y)i?(0()O)))+i?(0(o)i?(0(y)) 

-i?U(>O)i?U(O<)i?(0(>O)) 



So in the end we have the following counter term matrix 

/ 1 -i?(0(o)) -i?(</.(y)-<^(o)i?(0(o)) 



0- = 








1 




-R 



O)) 



1 




(61) 
(62) 
(63) 

\ 



-i?(0(y) -0(o<)i?(</'('O) 
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giving the counter terms for the graphs O, ^ and y)>i in its first row. The renormahzed character 
matrix 0+ follows from the matrix Birkhoff factorization (49) in Theorem 20, 0+ = According 
to the matrix- vector equation (50) of Theorem 20 we find 



{<t>V 



yielding the renormalized amplitudes for the graphs O, and S^k in its first row, which we 
write in transposed form 



0(0) -i? (0(0)) 

0( Y) - 0(o<)i?(0(o<)) - ^^(^l Y)) + r(^^{o)r{4^{o) 
-0(y)i2(0()O)) -0(o)i?(</'(y)) +0(o<)i?(0(o<)i?(</'(o< 
-R\ <j>[ f7\l -0(Y)i?(0()O)) -0(o<)i?(0(Y)) +0(o)i?(0()O)i?(</)(o))) 



Another example up to 3-loop order, including a graph with two disjoint 1-loop subdivergences 
is provided by taking X to be the set of graphs 



■^(4) := {ei := , 62 := |0<) , 63 := | Y)' ^4 := O)} 



(64) 



The new graph \^ is made of the two disjoint fish graphs O as subdivergences sitting inside of 
such a O graph. Remember that our graphs carry no external structure. The coproduct of this 
graph is given by 



+ 20 Y + oo ® <y. 



Let e Ga be the Feynman rules character. The coproduct matrix, respectively the character 
matrix are given by 



(f) :— (j) o Mjri 



(4) 



/I o Y £3 \ 

1 O o>o 

1 20 

V 1 / 



( 1 0(0) 0(Y) ^P) ^ 

1 0(0) 0(0)' 

1 20(0) 
\ 1 / 



In this example the only new counter term matrix entry we need to calculate is position (1,4) in 
0_. Applying formula (21) for its solution, we find 



/3i4 = -i?(^0(£l)^+i?(^0(o<)i?(0(O<)'))+it:(^0(Y)^(20(>O)) ) (65) 

-r{^{0) 2R(<l>{K>)R((t>{<>)))y 
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The identity i?(a)^ = 2R{aR{a)) — R{a'^) following from (1), and which is true only for commu- 
tative Rota-Baxter algebras, immediately implies 

Pi4 = -i?(^0(E3)) -i?(^0(O<)i?(<?>(>O))') +2i?(^0(Y)i?(0(O<))). (66) 

Likewise for the renormalized expression we find in entry (1,4) of the renormalized matrix 
respectively the 4th component of the vector 

((</.+ r)4 = 0(52)+0(>O)i?(</>(O<))'-20(Y)i?(0(>o)) (67) 

-rU(^) + 0(o<)it:(0(o<))' - 20( Y)it:(0()o))) . 



3.5 More examples and comments on matrix factorization 



As another illustration, let us consider the case of the truncated space 7(6) in (39) of undecorated 

rooted trees. For a given regularized character (f) : Tiq > A, we have the corresponding matrix 

in Equation (45) which we record below for easy reference. 







o M 



n 



(68) 



/ 1 





\o 



0(62) 0(63) 



1 








(e2j 
1 






0(64) 

0(63) 

0(62) 
1 






0(65) 

0(e2)^ 

20(62) 


1 




0(e6) \ 
0(e4) 
0(63) 
0(e2) 


1 / 



where Mt-c is the coproduct matrix in (41). Recall that the unit diagonal upper triangular matrix 
/3 0: = exp (-7^(x(^^))) = 1 + 7^ (exp*- (-x(^<^)) - l) (69) 
is solution to the equation (47) 

/3 = l-7^((0-l)/3). 

Remember that the second equahty in (69) follows from Proposition 32 respectively its matrix 
representation, with % in place of x, see Theorem 20. . 
Equation (47) is solved by the formula (21), and we find 



/ 1 -i?(0(e2)) 



(5:-- 










1 








/3l3 

-i?(0(e2)) 
1 






/3l4 
I32A 

-i?(0(e2)) 

1 






/3l5 
/?25 



1 





/9l6 
/526 
/936 

1 
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where, in abbreviating 0(ej) by e,, 1 < i < 6, 

Pi3 = -Rie;)+R{e2R(e2)), 

Pi4 = -Rie,)+R{e2Rie^))+R{esRie2))-R{e2Rie2Rie2))), 
/3i5 = -R{er,) + R{e2R{e2e2)) + R{e^R{2e2)) - R{e2R{e2R{2e2))) , 
/3i6 = -R(ee) + R{e2R{ei)) + R{e^R{e^)) + R{eiR{e2)) 

-R{e2R{e2R{es))) - R{e2R{esR{e2))) - R{e3R{e2R{e2))) + R{e2R{e2R{e2R{e2)))), 
P24 = -i?(e3) + i?(e2i?(e2)), 
/925 = -Riel) + R{e2R{2e2)) = R{e2)R{e2), 
P26 = -R{ei) + R{e2R{e^)) + R[e^R{e2)) - R[e2R{e2R{e2))) , 
(335 = -2i?(e2), 
PsG = -7?(e3) + i?(e2i?(e2)), 
= --^(62), and = 0. 

The counter term expressions for the graphs e^, i = 2, . . . , 6 are given in the first row of (3, 

(</._! = (lr|/9= (l,-i?(e2),-i?(e3) + i?(e2i?(e2)),/?l4,/9l5,/9l6)■ 
These calculations for the Pij match the results of Bogoliubov's counter term recursion applied to 
the coproduct matrix. For the renormalized matrix character = 0+ we find 

0+ = exp(^(x(^^))) = l-^(exp*^(-x(^^))-l) 
Then the renormalized expressions 0+(ej) for i = 2, . . . , 6 are obtained as components of the vector 

(0+|-(lr|0+=(lr|00- = (0|0- 

As a remark for the practitioner we mentioned that from Proposition 7 we derive the more familiar 
equation for the renormalized character 

$+ = 1 + ^((J-1)0_), (70) 

which is just Bogoliubov's classical R-operation giving the renormalized Feynman rules (matrix). 

3.5.1 Exponential approach of matrix factorizations 

In this subsection we leave the realm of combinatorial Hopf algebras of renormalization and come 
back to the results of Section 2.3, where we discussed the decomposition of upper triangular matrices 
with entries in a commutative Rota-Baxter algebra. 

The following dwells on the exponential approach to the calculation of the factor matrices a± 
in the factorization of a in Corollary 10. Other than its theoretical significance, it also relates to 
the exponential approach of the Birkhoff decomposition of Connes and Kreimer. 
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The reader should remember Theorem 9 asserting that with {A, R) being a commutative Rota- 
Baxter algebra, the triple {Ain{A),TZ,{Ain{A)k}k>i) forms a complete filtered Rota-Baxter alge- 
bra. 

Let us start with some properties of the exponential and logarithm functions for complete filtered 
algebras of upper triangular matrices defined in (13) respectively (14). For any 1 < < oo, a natural 
basis of M'^{A) is given by the matrices E^j e M'^{A), < i < j < n, where the entry at position 
(i, j) is 1, the rest zero. These matrices multiply according to EijE^i = Sj^En. We can express the 
logarithm and exponential in terms of this basis {Eij}. For a e d}ln{A), we have log(Q;) e Mn{A)i. 
So we have 

Za := log{a) = (aij) = ^ aijEy. 

0<i<j<n 

These a^j G A are called matrix normal coordinates (of the second kind). The concept of normal 
coordinates in the context of Connes-Kreimer renormalization theory appeared in [21]. 
For example, let the 3x3 matrix a in 9Jt3(A) be 

3 

a = ^ -Eii + aEi2 + bEi^ + c^as 
1=1 

with a,b,c G A. Note that a — 1 is strictly upper triangular and so (a — 1)'' = for k > 3. 
Therefore, we have 

Z^ := log{a)^a-l-^{a-iy^aEi2+(^b-^ac^Eis + cE23eMUA)i, (71) 
giving normal coordinates 

ai2 = a, q;i3 = 6 - -ac, a23 = c. 

Thus 

a = exp (^aEi2 + (^b - ^ac^ E13 + 0^23^ . (72) 

In general, for given a = {aij) G ^))tn{A) these matrix normal coordinates can be calculated by 
the formula 

aij = aij + ^ ^ -j-j-jau^ai^i^. . .ai^j, < i < j < n. (73) 

k=l i<h<l2<-<lk<j 

These new coordinates allow us to write any n x n-matrix a G 9Jl„(yl) as 



a = exp (Z„) = exp I ^ 

\o<i<j<n 



O-ij E^j 



In order to obtain the factorization a = Q!+Q!_^, we need the BCH-recursion (10) 

X{Z^) := Z^ - BCH (n{x{z^)),n{x{Za))) , 
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which allows us to calculate TZ{x{Za)) and iZ{x{Za)) in Ai'^{A)i. This is valid for any n < oo. 

We continue with our example of Ai'^{A). Note that Ai'^{A)^ — ior k > 3 and that commuta- 
tors of higher order in x are identically zero in Ai'^{A)i, due to the decreasing filtration. Thus by 
(17), we find 



= exp (-7^(x(^a)) 

= exp(^-7^(z,-^[:^(z,),7^(z,)])^ 

7^(z, - ^[^(z,),7^(z,)]) + i7^(z,)7^(z,) 



= 1 



We would like to underline, that up to this point we have not used the condition that 7?. is a 
Rota-Baxter map. Actually, for the factorization, we only needed that TZ + TZ — idM^{A), which is 
true for any linear map, and the BCH-recursion x (or x) ^ suitable topology, such as a complete 
filtered algebra. The Rota-Baxter structure only enters in the next step, when we replace the last 
term TZ{Za)TZ{Za) via Rota-Baxter relation. Keeping in mind that, for the Lie bracket [— , — ], we 
have 

[n{Za),n{z^)] = [Za - n{Za),n{Za)] = [Za,n{Za)], 



thus 



a_ = l-7^(z„) + -(7^(z,7^(z„))-7^(7^(z,)z,)) + 



= 1 

In matrix form we therefore get 



+^{n{n{z^)z^) + n{z^n{z^))-n{z^z^)) 
n{z^) + n{zMZa)) - InZaZa). 



(74) 
(75) 



/ 1 -R{a) ( - R{b) + lR{ac)) \ 



a. 



1 

\ 

/ 1 -R{a) 

1 

V 



/ R{aR{c)) \ 



+ 



I 



-R{c) 
1 

( - Rib) + R{aR{c))) \ 
-R(c) 

1 / 




\ 








/ lR{ac) \ 




\ 








(76) 



The inverse of «_ can be calculated, using the recursive formula (12) or directly from (17) 

= exp (7^(x(Z„)) 

= i + n(^z^- ^[n{z^),'R{z^)]^ + ^n{Za)n{Za) 



a_ 



l + 7^(z,) + 7^(7^(z„)z,)--7^( 

ZaZo 
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' 1 


R{a) 


R{h) - \R{ac) \ 




/ 


R{R{a)c 







1 


R{c) 


+ 










.0 





1 ) 




\ 







/ 1 


R{a) 


R{h) + R{R{a)c) - R{ac)) \ 









1 


R{c) 










Vo 





1 




) 






/ 1 


R[a) 


Rib) - R{R{a)c)) \ 











1 


R{c) 















1 


/ 







/ \R{ac) 


y 



We similarly calculate a+ = exp ^7?.(x('^a)) j and reach the following factorization for example 
of 3 X 3 matrix a e ^1{A). 



a — OL^OL 



-1 



/ 1 R{a) R{b) - R{aR{c)) \ 

1 R(c) 
Vo 1 ) 

This recovers the factorization in (20). 



f 1 R{a) R{b) - R{R{a)c) \ 
1 R{c) 

VO 1 ; 



(77) 



Wc finally make a remark on the normal coordinates for the example in (68). For the matrix 
representation of the character 0, ^a[0] = = exp (^Z^, with = log(cD), the strictly upper 

triangular matrix Z^ G AiQ(A)i follows by using the formula for the matrix normal coordinates 
(73) 



0<i<j<6 

^12-^12 + ^13-^13 + ^14-^14 + 015-E^15 + 016-£^16 + 
023-E'23 + 024-E'24 + 025-E'25 + 026-E'26 + 

034-E'34 + 035-E'35 + 036-E'36 + 045-E'45 + 046-£'46 + ^Se-E^SG 



(78) 



/o 


0(62) 


m) 


0(64) 


0(65) 


0(66) 


\ 








0(e2) 







0(64) 













m) 


20(^2) 


0(63) 



















0(62) 
























Vo 

















/ 



(79) 



We used the fact that G Ga is a character, i.e. an algebra homomorphism. We have the following 
simple polynomial expressions for e,, i = 1, . . . , 6, following from (73) 

1 



62 



62, 63 



63 - 2^262 



(80) 
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64 = 64 - 6263 + -626262, 65 = 65 - 6263 + -626262 (81) 
11 , , 

ee = 66 - 6264 - 2^363 + 626263 - -62626262. (82) 

These are exactly the rooted tree normal coordinates as they appear in [21], and [26, 29]. They can 
be calculated as well using the convolution product and the logarithmic map, = log* (id7^^)(6j), 
i = 1,...,6. 

We hope that these examples provided some insight into the underlying structure and cal- 
culational simplicity of the matrix factorization in the context of complete filtered Rota-Baxter 
algebras. 

Let us briefly summarize what we have found in this section. Upper triangular nxn matrices, for 
any 1 < n < oo, with entries in a commutative Rota-Baxter algebra {A, R), Ai'^{A) form a complete 
filtered Rota-Baxter algebra {Mn{A),TZ, {Ain{A)}k>i)- The complete filtration allows us to define 
a Baker-Campbell-Hausdorff based recursion relation, denoted by x '■ •^n(^)i ~^ •^n(^)i) which 
in turn gives rise to a decomposition of the group dJln{A) of upper triangular matrices with unit 
diagonal. The linear map TZ appearing in the definition of the recursion for x can be any linear 
map in End(A^^(A)). Choosing it to be a Rota-Baxter map gives rise to solutions of the matrix 
group factorization in terms of recursion equations. Theorem 20. 

3.6 Berg-Cartier's ansatz using the grafting operation on rooted trees 

Berg-Cartier [4] used a different but related approach to encode the derivations in in terms of 
lower triangular matrices. For this they made use of the pre-Lie insertion operation on Feynman 
graphs. 

For the sake of convenience, let us state briefly the definition of pre-Lie algebra. Let ^4 be a 
not necessarily associative K-algebra. Wc denote the multiplication : A ^ A ^ A in A hj 
concatenation, m^(a(8) 6) — a b, a,b & A. The associator is defined as (— , — , —)a : A x A x A ^ A, 

(a, b, c)a '■— a {b c) — (a b) c (83) 

for a,b,c E A. For A being an associative K-algebra we have (a, b, c)a = for all a,b,c & A. 

A (left) pre-Lie K-algebra (-P, o) is a K-vector space P, together with a bilinear pre-Lie product 
o : P X P — > P, fulfilling the (left) pre-Lie relation 

(a, 6, c)p = {b, a, c) p, Va, b,ce P. (84) 

or explicitly 

a o {b o c) — {a o b) o c — b o {a o c) — {b o a) o c, Va, b,c & P. 

The pre-Lie property is weaker than associativity, i.e. every associative K-algebra is evidently pre- 
Lie. The commutator [a, 6] a o b — b o a for a,b e P fulfills the Jacobi identity, making the 
K-vector space underlying P a Lie algebra. 
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In the rooted tree setting the process of insertion of Feynman graphs into other graphs becomes 
a grafting operation. The derivations (29) Zt form a pre-Lie algebra Zt' o Zf := '^^^^n{t',t"]t)Zt 
[16, 17, 19], which by anti-symmetrization defines the commutator of the Lie algebra g of derivations 
(30). This pre-Lie composition is used to define an action of the Zf on the vector space T freely 
spanned by the rooted trees (or Feynman graphs). The operator representing the action is denoted 
by s(t) : T —>■ T, for all t G T and defined as follows 

s,\t") ■.= J2n{t',t"-t) \t). 

As an example we calculate 

s(.)|.) = |I>, s(.)|I> = |i) + 2|A>, s(I)|.)= |J^, andby definition s(t) I Ir) := \t) \/ter. 

We used the ket-notation for the rooted tree vectors introduced in an earlier section. The rule for 
calculating the vector s{t) \t') G T is to graft the tree t in all possible ways to the tree t\ and to 
multiply each tree in the resulting linear combination by its symmetry factor. This action can be 
used to define a representation of the Z^'s in terms of lower triangular matrices, which are just 
the transposed of our upper ones (42). The difference between lower and upper triangular matrix 
representation reflects the fact that the former increases the degree by grafting trees, whereas the 
latter reduces them by "elimination" of subtrees. The matrix representation approach using the 
pre-Lie structure on Feynman graphs or rooted trees appears to be limited to representations of 
infinitesimal characters respectively characters. Relation (36) instead works for arbitrary elements 
in Hom(7^r,^)- 

It also appears that the ansatz chosen in [4] for the counter term matrix, denoted by Ci/^ in [4], 

(Ci/e)^' = exp ( - 5^ Cm)Ht)\ , (85) 

using the operator s{t) for t G T is not sufficient for several reasons. First, the representation 
of characters G Hom(?-^r, ^) via the exponential map demands the use of normal coordinates. 
Second, the counter term character 0_, or its matrix representation 0_, follows from the factor- 
ization of characters in the sense of Atkinson, see Theorem (4) respectively Spitzer's identity for 
non-commutative associative Rota-Baxter algebras. Theorem (3). Therefore one must include the 
particular properties of the Rota-Baxter relation as well as the BCH- recursion (7). This shortcom- 
ing in [4] becomes particularly evident when comparing Equations (18), (19) and the one following 
(19), therein (see Equation (67) above for the correct expression). Equation (18) is plagued with 
unwanted coefficients. The derivation of the expression after Equation (19) in [4] is problematic 
as it seems to assume the subtraction scheme map, denoted by C in [4], to be an idempotent al- 
gebra homomorphism. This is not true in general, e.g. in the MS scheme C = Rms, which keeps 
only the pole part of Laurent series, is a projector of Rota-Baxter type. Berg and Cartier ear- 
lier in their paper (page 18 in[4]) proposed a special rule for products in the image of the map 
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C, e.g. C{a)C{a) is supposed to read as \C{C{a)a), which in general is not sufficient to resolve 
the aforementioned inconsistency. Instead, in a commutative Rota-Baxter algebra (^4, Ft) we find 
\R{af = R{R{a)a) - \R{a^), a e A. 

The simple calculation of the counter term matrix (85) used in [4] would only apply to the 
subset of ladder trees (or Feynman graphs) denoted by Id C T, if normal coordinates were properly 
used^. This follows from the fact that linearly ordered ladder trees form a cocommutative Hopf 
subalgebra Ti^d C Tir, or equivalently, the dually defined derivations Zt^^ ~: Zn, indexed by the 
number of vertices of a ladder graph tgn form a commutative Lie subalgebra g^i C g, on which the 
BCH-recursion reduces to the identity map, xlgu — id- Also, a consistent motivation for reducing 
the algebraic Birkhoff decomposition of the character matrix (f) = (f)^ (pz} into a linear matrix- vector 
equation of the form (50) was not given. 

4 Conclusion 

We have shown that the Connes-Kreimer Hopf algebra approach to rcnormalization in pQFT can 
be entirely represented as a simple and efficient triangular matrix calculus. Decomposing an n x n 
upper triangular unital matrix as described above, provides us with the counter term matrix as well 
as the renormalized matrix. The matrix calculus allows for an efficient calculation of counter terms, 
and henceforth renormalized Feynman amplitudes. 
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